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Abstract 


This  paper  describes  a  new  procedure  for  estimating  parameters  of 
a  stochastic  activity  network  of  N  arcs.  The  parameters  include  the 
probability  that  path  m  is  the  longest  path,  the  probability  that 
path  m  is  the  shortest  path,  the  probability  that  arc  i  is  on  the 
longest  path  and  the  probability  that  arc  i  is  on  the  shortest  path.  V 
The  proposed  procedure  uses  quasirandom  points  together  with  information 
on  a  cutset  H  of  the  network  to  produce  an  upper  bound  of  0((log  K)N" 'H'+1/K) 
on  the  absolute  error  of  approximation  where  K  denotes  the  number  of 
replications.  This  is  a  deterministic  bound  and  is  more  favorable  than 
the  convergence  rate  of  1/K^  that  one  obtains  for  the  standard  error 
for  K  independent  replications  using  random  sampling.  It  is  also  shown 
how  series  reduction  can  improve  the  convergence  rate  by  reducing  the 
exponent  on  log  K  .  The  technique  is  illustrated  using  a  Monte  Carlo 
sampling  experiment  for  a  network  of  16  relevant  arcs  with  a  cutset 
of  H=7  arcs.  The  illustration  shows  the  superior  performance  of  using 
quasirandom  points  with  a  cutset  (plan  A)  and  the  even  better  performance 
of  using  quasirandom  points  with  the  cutset  together  with  series  reduction 
(plan  B)  with  regard  to  mean-square  error.  However,  it  also  shows  that 
computation  time  considerations  favor  plan  A  when  K  is  small  and  clan  B 
when  K  is  large. 

KEY  WORDS 
Longest  path 
Monte  Carlo  methods 
Network  analysis 
Quasi  random  points 
Shortest  path 


Introduction 


Although  comprehensive  procedures  have  been  developed  for  analyzing 
an  activity  network  with  deterministic  arc  completion  times  (e.g., 

Elmaghraby  1977),  no  comparable  set  of  comprehensive  procedures  exists 
for  a  network  in  which  the  arc  completion  times  are  random.  In  particular, 
a  stochastic  activity  network  is  properly  characterized  by  distributional 
information.  For  example,  on  independent  passes  through  such  a  network 
from  source  to  sink,  the  network  completion  times  on  each  pass  usually 
differ,  as  do  the  paths  identified  as  the  critical  paths.  Here,  the 
distribution  function  (d.f. )  of  network  completion  time  and  the  frequency 
with  which  each  path  is  the  critical  path  provide  useful  characterizations 
for  describing  network  behavior. 

Whereas  conventional  search  algorithms  have  been  the  principal  tools  of 
analysis  for  the  deterministic  case,  Monte  Carlo  methods  have  been  the 
principal  tool  for  the  stochastic  case.  For  example,  see  Elmaghraby  (1977) 
and  Sullivan,  Hayya  and  Schaul  (1982).  However,  the  performance  of  the 
Monte  Carlo  method  has  not  been  exceptional  until  recently  when  Sigal, 
Pritsker  and  Solberg  (1979)  showed  how  one  could  use  information  on  a 
cutset  of  a  network  to  improve  the  accuracy  of  estimate  for  a  given  cost. 
More  recently,  Fishman  (1983)  showed  how  the  cutset  information  can  be 
used  together  with  quasi  random  points  and  series  reduction  to  estimate 
the  d.fs.  of  network  completion  time  and  shortest  path  time  with  upper 

bounds  on  the  errors  of  estimate  that  converge  more  rapidly  than 

1  /2 

0(1/K  )  ,  the  rate  for  the  standard  error  in  a  Monte  Carlo  experiment 

using  random  sampling  on  K  independent  replications. 

The  purpose  of  this  paper  is  to  extend  this  combined  cutset/quasi  random 
points/series  reduction  approach  to  the  estimation  of  the  probability  that 
path  m  is  the  longest  path,  the  probability  that  path  m  is  the  shortest 
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I 


i 
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path,  the  probability  that  arc  i  is  on  the  longest  path  and  the 
probability  that  arc  i  is  on  the  shortest  path.  Section  1  introduces 
the  necessary  notation.  Section  2  describes  the  use  of  cutset  information, 
as  proposed  by  Sigal,  Pritsker  and  Solberg  (1979),  to  estimate  a  parameter 
by  Monte  Carlo  methods  with  less  variance  than  would  obtain  for  the  same 
number  of  replications  when  using  crude  Monte  Carlo  techniques.  Section  3 
then  describes  how  one  can  use  quasi  random  points  with  the  cutset  information 
on  each  replication  to  make  the  error  of  estimate  converge  more  rapidly  than 
in  the  case  of  random  sampling.  Section  4  extends  the  analysis  to  the 
remaining  estimators  of  interest.  Section  5  describes  how  series  reduction 
in  the  network  can  work  to  advantage  by  reducing  the  number  of  dimensions 
sampled.  Section  6  then  describes  an  experimental  design  for  evaluating 
the  benefits  to  be  derived  from  using  cutset  information  together  with 
quasi  random  points  and  series  reduction.  Section  7  presents  an  example 
which  assesses  the  benefits  of  each  sampling  plan  versus  its  relative 
computation  time. 

1 .  Definitions 

Consider  an  acyclic  directed  network  with  a  single  source,  single  sink, 

N  arcs  and  L  paths.  Let  X.|,...,XN  ,  the  passage  times  for  arcs  1,...,N  , 
be  independent  random  variables  where  has  distribution  function  (d.f.) 

Fi  on  [0 ,°°)  and  inverse  distribution  function  G^u)  =  min[x:  F^x)  >  u  , 

0  <  u  <  1]  .  The  completion  time  of  path  m  is 


T  =  r  a.  X.  =  V.  _  X. 

m  “ • im  i  Mel  i 

i-l  m 


(1) 


where  a.  =1  if  arc  i  is  on  path  m  ,  a.  =0  otherwise  and  I  denotes 
im  im  m 

the  set  of  arcs  on  path  m  .  Consider  the  network  completion  time 


T  =  max{T.| , . . .  ,T^) 


(2) 


and  the  shortest  path  time 


T*  =  min(T-| ,. . .  J^)  .  (3) 

Characterizations  of  the  network  include: 


a. 

pr(T  *  t) 

0  *  t  <  » 

b. 

★ 

E  T 

c. 

gm=pr(m  is  the  longest  path) 

m=l ,. . .  ,L 

d. 

r^=pr  (arc  i  is  on  the  longest  path) 

i=l , . . .  ,N 

e. 

pr(T*  <;  t) 

O 

IA 

r+- 

A 

8 

f. 

E  T* 

g. 

hm=pr(m  is  the  shortest  path) 

m=l , . . .  ,L 

h. 

s^=pr(arc  i  is  on  the  shortest  path) 

i=l , . . .  ,N 

Fishman  (1983)  describes  a  highly  efficient  method  of  estimating  a,  b,  e 
and  f  based  on  cutset  information,  quasi  random  points,  and  series 
reduction.  The  present  paper  extends  the  method  to  the  estimation  of  the 
probabilities  {gm  ;  m=l,...,L)  ,  {hm;  m=l,...,L}  ,  {ri  ;  i=l,...,N} 
and  {s..;  i=l,...,N>  . 

2.  Cutsets 

Observe  that 

gm  =  pr[Tm  >  Tn;  njhn,  n=l,...,L]  m=l,...,L  .  (4) 

Because  of  the  shared  arcs  among  paths,  T^,...,TL  are  usually  not 
independent  and  therefore  gm  generally  has  no  readily  accessible  analytical 
form.  However,  one  can  make  use  of  the  cutset  information  as  in  Sigal, 
Pritsker  and  Solberg  (1979)  to  derive  a  useful  expression  for  estimating  g^  . 

Let  H  denote  a  cutset  of  the  network.  As  used  here,  a  cutset  denotes 
a  set  of  arcs  that  connects  a  set  of  nodes  W  containing  the  source  with 


a  set  of  nodes  W  containing  the  sink.  Also  assume  that  each  path  has 
only  one  arc  in  H  .  If  each  arc  in  H  points  from  W  to  W  ,  H 
is  called  a  uniformly  directed  cutset  (UDC).  We  assume  that  H  is  a  UDC. 


Let 


and 


Y. 


l 


sup 
m=l , . 


m  1 1  •  ( .  f  L 


ieH  (5) 


where  i  =1  D H  denotes  the  arc  on  path  m  that  is  in  the  cutset, 
mm 

Then  X1+Yi  is  the  maximal  time  of  all  paths  that  use  arc  i  ieH  and 
pr(m  is  the  longest  path | ,. . . .S^) 


Cnu„  F^t-YpUF^lt-S.) 

W. 


m 


m 


=  0 


Therefore, 


if 


otherwise. 


(6) 


m 


n 


ieH 

i7i 


Fi<VYi> 


if 


T  = Y.  +X. 
m  i 

m  m 


m 


(7) 


=  0  otherwise 

has  expectation  gm  .  Let  the  subscript  j  denote  replication  j  .  Then 
for  K  replications 


where 


rfiK  "  K 


1  vK 
K  L.  . 
J=1 


6(1 


mj 


v' ,n1;«  F<(W 

lj*l 

m 


(8) 


<5  ( X )  =  1 
=  0 


if  x=0 


otherwise. 
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is  an  unbiased  estimator  of  gm  .  If  the  replications  are  statistically 
independent,  then  var  is  proportional  to  1/K  . 

In  principle,  (8)  is  the  estimator  of  gm  that  Si  gal ,  Pritsker 
and  Sol  berg  propose.  They  note  that  it  should  give  smaller  variance  for 
K  independent  replications  than  the  more  conventional  estimator 


v  =  «<YV 


where 


Tj  ■  . V 


However,  both  var  §m(<  and  var  gmj<  converge  as  l/K  .  In  this 
paper  we  propose  the  use  of  quasi  random  points  to  speed  up  the 
convergence  of  the  error  of  estimate.  See  Niederreiter  (1978)  for  a 
comprehensive  review  of  the  associated  theory. 

3.  Quasi  random  Points 

Suppose  one  wants  to  estimate  the  multivariable  integral 

0  =  j  ...j  f (x, , . . .  >xm)  dx-j...dx^,  (9) 

0  0 

by  the  summation 

9K  =  1  £j=l  f(ulj”'-,UNj) 

where  {u^,. . .  ,u^. :  j=l,...,K)  is  a  sequence  of  points  uniformly 

distributee!  in  the  N-dimensional  unit  hypercube  Q^=[0,1]  .  Let  A(R;K) 

denote  the  number  ot  joints  Uj=(u^. ,. . .  »u^j)  j=l,...,K  that  fall  in  the 

region  R  C  •  Then  u^,...,uK  are  said  to  be  uniformly  distributed 
in  if 

, .  A(R;K)  n 

hyM  - 1 

for  all  R={(x] ,. . . ,xN) :  aj  s  x.  <  3.  ;  i=l,...,N}  .  Here  x(R)  is  the 
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measure  of  volume  of  R  in  .  We  define  the  extreme  di screpancy 
of  u=ur...,uK  as 

D.  =  DJN.u)  =  sup  |MMI  _  X(R)  (10) 

K  K  K 

for  R={  ( x-j , . . .  ,x^) :  0  ^  x^  <  8^  ;  i=l,...,N}  .  If  f  is  of  bounded 

variation  in  the  sense  of  Hardy  and  Krause,  then 

|eK-e|<  V(f,N)DK  .  (11) 

where  V(f,N)  is  a  function  of  the  bounded  variation  of  f  in  N  and 
lower  dimensions.  Most  importantly,  uniform  sequences  exist  for  which 
Dk  <  0( (log  K)N/K)  .  (12) 

This  is  an  impressive  result,  for  in  the  special  case  in  which 

{u . . :  i=l,...,N;  j=l,...,K}  is  a  sequence  of  i.i.d.  random  variables 

*  J 

uniformly  distributed  on  [0,1)  ,  one  has  only 

J2K  D 

lim  sup  — =  1  w.p.  1,  (13) 

K-h»  4 1  og  log  K 

a  result  due  to  Chung  (1949)  for  N=1  and  to  Kiefer  (1961)  for  N>1 
Hal  ton  (1960)  describes  a  uniformly  distributed  sequence  for  which 
(12)  holds.  If  R  >  2  is  an  integer,  then  every  non-negative  integer  n 
has  an  expansion  of  the  form 

n  =  I^q  a-j^  a^e(0,l , . . .  ,R-1 )  p^.) 

0  <  i  <  m  and  m=LlogRnJ  . 

Corresponding  to  this  expansion  one  has  the  radical  i nverse  function 

$R(n)  =  I^q  a.R-1"1  .  (15) 

Then  Halton  shows  that  the  sequence  Ur,  ( j <f>D  (j);  j=0,l ,. . . ,K-1 }  , 

K1  KN 

where  R^,...,R^  are  pairwise  relatively  prime,  gives  the  bound  (12). 
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We  now 

specialize  these  results  to 

the  network  analysis  at  hand. 

Define  for 

j=l , . . . ,K 

xio  =  Gi(uO> 

i=l , . . .  ,N 

(16) 

Vi  "  aircXi  j 

m=l , . . .  ,L 

(17) 

Y . .  =  sup  T  .  -  X.  . 

m~  i  ^  *  •  •  j  l  m 

ieH  . 

(18) 

Then  it  is  not  difficult  to  see  that  *  can  write  gm  as  a 

(N-  |H|+1  )-dimensional  integral  over  r  ,^|+-j  with  integrand  of  bounded 
variation  and  that  (8),  using  (16),  (  'd  (18),  is  an  approximation 
to  this  integral.  As  a  result,  one  has  tor  a  given  network  using  the 
Hal  ton  sequence  the  deterministic  bound 

IV-9„!  S0((lo9K)N-lHl+1/K)  ,  (19) 

which  makes  clear  the  desirability  of  choosing  as  H  the  cutset  of 

maximal  cardinality. 

4.  Estimating  {r . >  ,  {h  }  and  { s . } 

_ 3  i  m _  l 

This  section  describes  estimators  analogous  to  (8)  for 

r^  =  pr  (arc  i  is  on  the  longest  path)  i=l,...,N  (20) 

hm  =  pr  (m  is  the  shortest  path)  m=l,...,L  (21) 

and 

si  =  pr  (arc  i  is  on  the  shortest  path)  i=l,...,N.  (22) 

Before  doing  so,  we  explain  the  motivation  for  studying  {r^}  and  { s ^ }  . 

If,  say,  r.  turns  out  to  be  considerably  greater  than  the  remaining  arc 
probabilities  on  the  longest  path,  then  there  is  good  reason  to  believe 
that  the  completion  time  associated  with  arc  i  is  a  principal  determinant 


of  the  longest  path  time.  This  knowledge  of  the  relative  importance  of 
arc  i  may  lead  to  a  modification  of  the  activity  associated  with  this 
arc  for  the  purpose  of  shortening  the  longest  path  time.  Analogous 
considerations  apply  to  {s.}  .  Since 

r .  =  VL  ,  a .  g 
1  Lm=\  imam 

one  has  the  estimator 


A 

r . 


MT1=1  * 


iK  imymK 

Moreover,  it  is  easily  seen  that  |r.„  -  r. !  <  0(log  K)^' '^‘+1/K)  for 

1  f\  1 

the  Halton  sequence. 

The  estimation  of  {hi  follows  closely  that  for  (q  }  .  Let 

m 


Z.  =  inf  S 
i  m 

m=l , . . .  ,L 

a. 

lm 


i  eH 


so  that  X.+Zi  gives  the  minimal  time  of  all  paths  that  use  arc  i. 
Then 

pr(m  is  the  shortest  path  | .... ,SL) 

=  f{UUH  n-Fi(t-Zi)])dFi  (t-Sm) 


i/i 


if  Z.  =S 
l  n 


otherwi se . 


Again,  let  the  subscript  j  denote  replication  ,1  .  Also  let 


mj 


"  1  *  i  <  H 

i/i 


[1-F.(T  ,-Z.  .)] 
i  mj  ij  J 


if  T  ,-Z.  .  +  X . 


V 


'mJ 


=  0 


otherwi se . 
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so  that  our  estimator  has  the  form 

*,K  -  Cl  •  (28) 

with  error  |s.K-s.|  s  0((1og  K)N_|H|+1/K)  for  the  Halton  sequence. 

5.  Series  Reduction 

As  shown  in  Section  2,  the  rate  of  convergence  of  the  error  of 
approximation  (10)  is  bounded  by  a  function  0((log  K)N/K)  .  Using  the 
cutset  H  enables  one  to  bound  this  convergence  more  tightly  by 
0((log  K)^  |H|+1/K)  and,  presumably,  using  the  cutset  H  with  maximal 
cardinality  is  best.  In  fact,  any  preanalysis  that  leads  to  a  smaller 
exponent  for  log  K  is  desirable  from  the  viewpoint  of  a  tighter  error 
bound.  The  technique  of  series  reduction  is  one  such  approach  to 
further  reduction. 

Figure  1  shows  an  activity  network  due  to  Rattersby  (1970)  for  the 
partial  overhaul  of  a  unit  in  an  oil  refinery.  Table  1  shows  the  corresponding 
incidence  matrix.  Note  that  arcs  17  and  18  are  both  on  all  eight  paths. 

Insert  Fig.  1  and  Table  1  about  here. 

Therefore,  we  can  ignore  these  arcs  and  reduce  the  number  of  relevant  arcs 
for  analysis  from  18  to  N= 1 6  arcs  without  any  effect  on  our  results. 
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Observe  that  e^,  eg,  eg,  eg,  e-j q  and  form  a  cutset  H  for  the 
netwonc.  Since  N=16  arcs  and  |H|=7  the  exponent  for  log  K  is 
N- | tf 1+1=10  ,  Now  note  that  arcs  and  occur  only  on  paths  1  and 
2  ,  eg  and  e^g  occur  only  on  path  5  ,  e^  and  e^  occur  only  on 
paths  6,  7  and  8  and  eg  and  e^  occur  only  on  paths  7  and  8  . 
Suppose  we  now  replace  these  arc  pairs  by  single  arcs  e,+e.  ,  e.+e,,-  , 
e7+e16  and  eg+e^  .  Then  the  network  has  N *  =1 2  arcs  and  the  cutset 
H'  with  e2  ,  eg  ,  eg  ,  eg+e^g  ♦  eg  ,  e^  and  e^  has  cardinality  7 
so  that  N'-|H' |+1=6  is  the  resulting  exponent  for  log  K  in  the  bound 

A  A 

on  l9m|<-9n^  anc*  •  Note  t*iat  if  one  adopts  this  series 

reduction  it  is  not  possible  to  estimate  r..  and  s^  for  those 
original  arcs  in  the  reduction. 

The  feasibility  of  implementing  this  series  reduction  in  practice 

depends  on  whether  or  not  the  d.fs.  of  the  summed  arc  times  are  convenient 

for  computation  and  on  the  relative  computation  costs  with  and  without 

series  reduction.  Suppose  Wj,...,W  are  independent  exponentially 

distributed  random  variables  with  distinct  rates  0-|,...,6n  .  Then  for 

n>l  V  =W.+...+W_  has  the  d.f. 
n  1  n 


fv  (O  =  Ii=l  ni  =  ](i  Bj/8i  )_1e 
n  W 
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(29) 


For  summed  arcs  in  the  cutset  one  need  only  replace  their  individual  d.fs. 
in  (8)  and  (26)  by  the  corresponding  d.f.  from  (29).  For  summed  arcs  not  in 
the  cutset  one  needs  to  solve  the  equation 


v  ■  fv  (n 

n 


(30) 


■i 
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for  T*  ,  v  being  the  appropriate  element  of  the  Halton  sequence.  This 
can  be  done  readily  by  the  Newton-Raphson  method.  A  suggested  starting 
point  for  the  iteration  is 

TJ  =  -t1n(l-*)]IJ.,  ^  .  (31) 

We  evaluate  this  approach  in  our  example  where  we  study  the  relative 
benefits  of  the  cutset  information,  quasi  random  points  and  series  reduction. 
As  we  show,  the  relative  benefit  of  series  reduction  becomes  more  apparent 
as  the  sample  size  K  increases.  In  fact,  for  small  K  series  reduction 

may  be  less  preferable. 

6 .  Experimental  Design 

This  section  describes  the  layout  for  an  experiment  designed  to 
determine  the  extent  to  which  quasi  random  points  lead  to  accelerated 
accuracy  when  estimating  {g^}  and  {h^}  for  the  network  in  Fig.  1. 

In  particular,  we  introduce  a  degree  of  randomness  into  the  experiment 
in  order  to  compute  estimates  of  the  mean-square  errors  of  our  point 
estimates.  We  then  study  the  behavior  of  these  sample  mean-square  errors 
as  K  increases.  Note  that  we  choose  this  method  of  evaluating  our  proposed 
method  merely  because  it  is  necessary  to  find  a  common  basis  for  comparing 
results  based  on  a  deterministic  method  with  baseline  results  using  random 
sampling. 

Consider  an  experiment  consisting  of  Q  statistically  independent 
blocks  or  macro repl ications  each  of  K  mi crorepli cations.  Let  (vim: 
i=l,...,N;  m=l,...,Q}  denote  a  sequence  of  i.i.d.  random  variables 
from  u[0,1)  and  define 

Cim  =  (c:  *R.(c)  =  vim)  i=1---*N  m=l,...,Q.  (32) 


Then  on  macroreplication  m  we  use  the  quasi  random  point 

+  J-1 ) »•  •  •  »4>r  (CNm  +  j-1 )  on  microreplication  j  for  j=l . K 

Let  0jm  denote  an  estimate  of  a  particular  quantity  e 
computed  on  microreplication  j  on  macroreplication  m  .  Then 
01 . ,6q  ,  where 

5m  =  1  &1  8jm  ra'1 . <>•  (33) 

are  i.i.d.  random  variables  with  sample  mean-square  error 

SK  =  (pf  Cl  (§m  '  eQ)2  (34) 


where 


eQ  =  Q  Cl  6m  * 

2 

It  is  the  behavior  of  sK  versus  K  that  interests  us.  Note  again  that 
the  only  reason  for  randomizing  the  seeds  (vim  ;  i=l,...,N}  on  each  of 
the  Q  macroreplications  is  to  enable  us  to  derive  the  estimates  of 
mean-square  error. 

To  provide  an  instructive  normalization  we  also  run  M  statistically 
independent  macroreplications  each  of  K=1  microreplication.  Let 

1  M 

denote  the  resulting  estimates  of  e  each  with  sample  mean-square  error 

WM  =  Cl  "  V2  (35) 

where 


0M  =  i  I  t  . 

M  M  cm= 1  m 


2  2 

Then  the  quantity  wM  /  Ks^  should  increase  as  K  increases  if  accelerated 
convergence  is  occurring. 
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7.  Exampl e 

To  illustrate  the  proposed  procedures  we  use  the  network  of  Fiq.  1 
for  which  we  estimate  (gm;  m-1,...,8}  and  {hm;  m=l,...,8}  .  Monte 
Carlo  experiments  were  performed  using  quasirandom  points  and  the 
cutset  H'={e2,  e3,  eg,  eg  +  elg,  eg,  e1Q,  e^}  with 
Plan  A 

No  series  reduction  on  the  remaining  arcs 

{el’  e4*  e7*  e8*  ell»  el 3 *  e14*  e16} 
and  then  with 

Plan  B 

Series  reduction  on  the  remaining  arcs 

{el  +  e4,  e?  +  e]6,  eg  +  e^,  e]3,  e]4  } 

Motivation  for  this  choice  of  plans  arose  from  the  observation  that 
plan  A  requires  (29)  to  be  evaluated  once  per  microreplication  whereas 
plan  B  additionally  requires  evaluation  of  (29)  and  its  derivative 
iteratively  (by  the  Newton-Raphson  method)  for  each  of  the  three 


reductions  e^  +  e4  ,  e^  +  e^g  and  eg  +  e^ 


While  it  is  true 


that  plan  B  also  requires  three  less  quasirandom  points  per  micro¬ 
replication,  the  relative  cost  for  an  iterative  computation  needs 
to  be  considered  explicitly. 

Expression  (8)  and  (26)  were  used  to  estimate  each  g^  and 
respectively.  For  the  baseline  case  of  pure  random  sampling  gm  and  h 
were  estimated  by 


1  rM  /  ^ 

gmM  =  M  Ej*l  6(Tj  '  V 


^mM  =  M  £j=l  6(T*j  "  W  ’ 


respectively,  where 


T*j  =  min^Tl j  * " ' ' *TLj ) 


Tj  =  max^Tl j ’ *  *  * ,TLj) 


Note  that  estimates  of  {r. }  and  {s^ }  could  in  principle  also  be 

computed  at  small  marginal  cost. 

2  2 

Let  w^(p)/K  and  s^p)  denote  the  sample  mean-square  errors  for 
the  estimators  of  p  .  Then  Table  2  lists  for  selected  K  the  averages 

Bh(W)  =  1 

for  the  shortest  path  using  sampling  plans  W=A  and  W=B 


Insert  Table  2  about  here. 


in  columns  3  and  4,  respectively,  and  the  averages 


6g(H)  "  5  ^1[wM(9m)/KsK(9m)1 
for  the  longest  path  in  columns  5  and  6  respectively. 

These  results  are  most  encouraging  for  they  indicate  substantive 
increases  in  these  ratios  as  the  number  of  microreplications  K  increases 
Also,  they  show  the  relative  desirability  of  using  series  reduction  as 
K  increases.  Note  that  computation  time  considerations  limited  the 
number  of  macroreplications  for  K  >  2^  . 

When  evaluating  a  proposed  accuracy  accelerating  technique  one  also 
needs  to  study  its  cost  relative  to  the  cost  of  crude  Monte  Carlo  sampling 
One  way  to  make  this  assessment  is  in  terms  of  the  CPU  times  required 
per  observation  with  each  sampling  plan.  To  simplify  the  discussion. 
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let  us  concentrate  on  the  estimation  of,  say,  gm  .  Suppose  that  a 
single  observation  using  pure  random  sampling  has  mean-square  error  V  . 
Then  over  K-j  independent  observations  or  microreplications  the  resulting 
estimate  of  gm  has  mean-square  error  VK  =V/K^  .  Suppose  that  for  K 
observations  of  the  quasi  random  point  approach  the  estimator  of  gm  has 
mean-square  error  .  Let  denote  the  ratio  of  cost  per 
microreplication  for  plan  W  using  quasirandom  points  over  the  cost 
per  microreplication  using  independent  random  sampling.  If  we  begin 
by  taking  K  observations  using  plan  A  and  observations  using 

independent  random  sampling,  then  each  experiment  takes  exactly  the 
same  amount  of  time.  Moreover,  the  ratio  of  their  mean-square  errors  is 


VK]  \I/K} 
fiK  fiK 


1  .  V/K 
“A  QK 


Let  $.  denote  the  value  of  the  estimate  on  independent  observation 

2  2 

i=l,...,M  .  To  estimate  V  we  use  wM  so  that  wM/K  is  the  estimate  of 
mean-square  error  V/K  for  the  case  of  pure  random  sampling  for  any  K  . 
Also  ls£  :  K=2^  -,  j=7,...,14l  are  the  estimates  of  mean-square  errors 

for  plan  A  for  the  selected  values  of  K  . 

To  estimate  the  average  ratios  of  mean-square  errors  for  equal  run 

times  (using  we  use  8^(W)/uw  and  8g(W)/w^  for  sampling  plans 

W=A,B  .  For  the  case  of  pure  random  sampling  one  microrepl i cation 
required  an  averaqe  of  7.37  microseconds  on  an  IBM  4341 /M2  whereas 
plan  A  required  an  average  of  12.25  microseconds  per  microreplication. 

Plan  B  required  an  average  of  14.85  microseconds  per  microreplication. 
Therefore,  l/u)^=7.37/12.25=.50  when  comparing  random  sampling  to 
plan  A  and  1/wg=7. 37/14.85=. 49  when  comparing  random  sampling  to 
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pl an  B.  Columns  7  through  10  show  eh(A)/u>A  ,  ( B )/ojb  ,  6  ( A)/oja 

and  6„(B)/u0  ,  respectively  for  selected  numbers  of  micro- 

9  o 

replications  K  . 

These  results  are  sobering  for  they  show  that  the  relative  desirability 
of  series  reduction  in  the  noncutset  arcs  does  not  set  in  until 
for  {hm>  and  for  (gm>  .  Let  us  now  put  this  cost  analysis 
into  perspective.  Since  and  oig  are  essentially  independent  of  K 
and  gene-ally  uig  >  ,  there  is  always  a  crossover  point  K*  for  each 

network  at  which  one  prefers  plan  A  to  plan  B  whenever  K  <  K*  and 
prefers  plan  B  to  plan  A  whenever  K  >  K*  .  Since  for  an  arbitrary  network 
one  does  not  know  K*  ,  the  choice  between  the  two  plans  is  not  trivial. 
Therefore,  the  advice  to  do  series  reduction  whenever  possible  cannot  be 
accepted  unqualifiedly  unless  the  planned  number  of  microreplications  K 
is  substantially  large. 
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